Noise invoked resonances near a homoclinic bifurcation 
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Stochastic Resonance (SR) and Coherence Resonance (CR) have been studied experimentally 
in the discharge plasma close to a homoclinic bifurcation. For the SR phenomena, it is observed 
that a superimposed subthreshold periodic signal can be recovered via stochastic modulations of the 
discharge voltage. Furthermore, it is realized that even in the absence of a subthreshold deterministic 
signal, the system dynamics can be recovered and optimized using noise. This effect is defined as 
CR in the literature. In the present experiments, induction of SR and CR are quantified using the 
Absolute Mean Difference (AMD) and Normalized Variance (NV) techniques respectively. AMD is 
a new statistical tool to quantify regularity in the stochastic resonance and is independent of lag. 

PACS numbers: 52.80.Hc 05.40.-a 05.45.Xt 05.20.-y 



CZ3 

11: 

O 
>v 

: 



> 

o 

cn 
o 

oo 
O 



X 



I. INTRODUCTION 

Stochastic Resonance (SR) is a phenomena in which 
the response of the nonlinear system to a weak periodic 
input signal is amplified/optimized by the presence of 
a particular level of noise [Ij, i.e, a previously untrace- 
able subthreshold signal applied to a nonlinear system, 
can be detected in the presence of noise. Furthermore, 
there exists an optimal level of noise for which the most 
efficient detection takes place [1, SR has been ob- 
served in many physical, chemical and biological systems 
S B 0, 0, M, El- Coherence Resonance (CR) 
is the phenomena wherein regularity of the dynamical 
behavior emerges by virtue of an interplay between the 
autonomous nonlinear dynamics and the superimposed 
stochastic fluctuations. In CR, analogous to SR, the ex- 
tent of provoked regularity depends upon the amplitude 
of added noise 0, U The CR effect too has been studied 
exhaustively, both theoretically and experimentally, in a 
wide range of nonlinear systems [jj, El EI EE El E3 • 

The nonlinearity in plasma systems arises from the 
most fundamental processes, namely the wave-wave and 
wave-particle interactions. Different modes may be ex- 
cited due to nonlinear coupling of waves and plasma com- 
ponents and the character of the oscillations is primar- 
ily determined by the plasma parameters and perturba- 
tions [ll| QjJ [2(|. In the present work, possibility of 
observing noise invoked resonances in a glow discharge 
plasma is explored. However, as a precursor to the ex- 
periments involving noise, a systematic analysis of the 
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autonomous dynamics is performed. This includes iden- 
tification and characterization of the bifurcation in the 
vicinity of the set-point employed for the noise related 
experiments. 



II. EXPERIMENTAL SETUP 

The experiments were performed in a hollow cathode 
dc glow discharge plasma. The schematic diagram of the 
experimental setup is presented in Fig[T] A hollow Stain- 
less Steel (S.S) tube of length « 7 cm and of diameter 
(<f>) w 45 mm was used as the cathode and a central rod 
of length w 7 cm and <f> w 1.6 mm was employed as the 
anode. The whole assembly was mounted inside a vac- 
uum chamber and was pumped down to a pressure of 
about 0.001 mbar using a rotary pump. The chamber 
was subsequently filled with the Argon gas up to a pre- 
determined value of neutral pressure by a needle valve. 
Finally a discharge was struck by a dc Discharge Voltage 
(DV), which could be varied in the range of 0—1000 V. 



The noise and subthreshold periodic square pulse gener- 
ators were coupled with DV through a capacitor [Fig. Q] . 
In all the experiments DV was used as the bifurcation 
parameter while the remaining system parameters like 
pressure etc., were maintained constant. The system ob- 
servable was the electrostatic floating potential, which 
was measured using a Langmuir probe of diameter <f> = 
0.5 mm and length 2 mm. The tip of this Langmuir 
probe was placed in the center of the electrode system 
as indicated in Fig. [TJ The plasma density and the elec- 
tron temperature were determined to be of the order of 
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FIG. 1: Schematic diagram of the cylindrical electrode sys- 
tem of the glow discharge plasma. The probe was placed at 
a distance I « 12.5 mm from the anode. Signal and noise 
sources were coupled to the Discharge Voltage (DV) through 
a capacitor. 



10 7 cm~ 3 and 3—4 eV respectively. Furthermore, the elec- 
tron plasma frequency (f pe ) was observed to be around 
28 MHz, whereas the ion plasma frequency (/pi) was mea- 
sured to be around 105 kHz. 



III. AUTONOMOUS DYNAMICS 

Before studying the noise induced dynamics, we char- 
acterized the behavior of the autonomous system. Not 
surprisingly, it was observed that at different chamber 
pressures, discharge struck at different voltages. Fig [2] 
shows the breakdown voltage (T4r) at different pd, where 
p and d are the filling pressure and radius of the cathode 
respectively. This breakdown voltage (Hr) initially de- 
creases with an increase in pd, goes through a minimum 
value resembling a typical Paschen curve and then begins 
to increase with increasing pd. It is observed that the sys- 
tem is excitable for the region pd > Paschen minima. In 
the lower side it shows Self Organized Criticality [2l[ . In 
this excitable domain, the system dynamics are irregular 
(complex) at the initial stages of the discharge voltage 
and upon increasing DV they become regular (period- 
one) as shown in Fig [3J Further augmentation of DV 
modifies the oscillation profile and results in the induc- 
tion of typical relaxation oscillations {22} ■ 
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FIG. 2: Vbr vs pd (Paschen curve) for our experimental sys- 
tem. The minimum occurs at (1.69 mbar-mm, 251 V). The 
system is excitable for the pd > minimum of the curve. 



fixed point behavior. Time traces from top to bottom 
in the left panel of Fig [4] depict this period lengthening 
of the oscillatory behavior. A systematic analysis of the 
increment in the period (T), presented in Fig.|¥[(a) right 
panel] , indicates that the autonomous dynamics undergo 
a critical (exponential) slowing down. Consequently, the 
ln|V — Vh\ vs T curve can be fitted by a straight line, 
where Vjj is the bifurcation point separating the oscilla- 
tory domain and the steady state behavior. The results 
of Fig [5] indicate that the system dynamics undergo a 
homoclinic bifurcation at Vh resulting in the loss of os- 
cillations. 



An anode glow is observed with these oscillations. Figs 
EJa) shows that the glow with largest size, appears when 
the discharge is struck at a typical pressure of 0.95 mbar 
and its size decreases with increase in the DV until it 
finally disappears [Figs[5ja)-l5jh)]. This may some types 
of unstable structure in the plas ma and produces such 
oscillation of the instabilities (231 . 



The time period (T) of these relaxation oscillations in- 
creases dramatically upon further incrementing DV. This 
eventually results in the vanishing of the limit cycle be- 
havior beyond a critical DV (Vh)- For larger values of 
DV, the autonomous dynamics exhibit a steady state 



IV. NOISE INDUCED RESONANCE IN GLOW 
DISCHARGE 

In this section experimental results involving noise gen- 
erated resonances, namely SR and CR are presented. 
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FIG. 3: Sequential change in the raw signal (normalized) at 
0.89 mbar for different voltages: (a) 288 V; (b) 291 V; (c) 295 
V; (d) 301 V; (e) 304 V; (f) 307 V; (g) 327 V; (h) 385 V; (i) 
466 V. All y-axes range form -1 to 1. 



A. Stochastic Resonance 

For our experiments on stochastic resonance, the refer- 
ence voltage Vo was chosen such that Vo > Vh and there- 
fore the autonomous dynamics, by virtue of an underly- 
ing homoclinic bifurcation, exhibit steady state behav- 
ior. The discharge voltage V was thereafter perturbed 
V = Vo + S(t) + D£, where S(t) is the subthreshold pe- 
riodic pulse train chosen for which V = Vo + S(t) > Vh, 
(subthreshold signal does not cause the system to cross 
over to the oscillatory regime) and D£ is the added Gaus- 
sian white noise £ with amplitude D. Subthreshold pe- 
riodic square pulse of width 20 (is and duration 2 ms 
was constructed using Fluke PM5138A function genera- 
tor. Meanwhile, the gaussian noise produced using the 
HP 33120A noise generator was subsequently amplified 



using a noise amplifier. 

Fig. [(Ja) — (c) show time series of the system response 
in the presence of an identical subthreshold signal for 
three different amplitudes of imposed noise. The sub- 
threshold periodic pulse train is also plotted, in the top 
most graph of the left panel, for comparison purposes. 
Fig. [UJa) shows that there is little correspondence be- 
tween the subthreshold signal and the system response 
for a low noise amplitude. However, there is excel- 
lent correspondence at an intermediate noise amplitude 
[Fig. [6][b)]. Finally at higher amplitudes of noise the 
subthreshold signal is lost amidst stochastic fluctuations 
of the system response [Fig. |5tc)]. Absolute mean differ- 
ence (AMD), used to quantify the information transfer 
between the subthreshold signal and the system response, 
is defined as AMD = abs(mean(-g- — 1)). t v and 8 are the 
inter-peak interval of the response signal and mean peak 
interval of the subthreshold periodic signal respectively. 

Fig [6](d) shows that the experimentally computed 
AMD versus noise amplitude D curve has a unimodal 
structure typical for the SR phenomena. The minima 
in this curve corresponds to the optimal noise level for 
which maximum information transfer between the input 
and the output takes place. 



B. Coherence Resonance 

For the experiments on coherence resonance DV (Vo) 
was located such that the floating potential fluctua- 
tions exhibit fixed point behavior. In order to min- 
imize the effect of parameter drift, a set-point (Vq) 
quite far from the homoclinic bifurcation (Vh) was cho- 
sen. Subsequently superimposed noise on the discharge 
voltage was increased and the provoked dynamics ana- 
lyzed. The normalized variance (NV) was used to quan- 
tify the extent of induced regularity. It is defined as 
NV = std(tp) I mean(t p ) , where t p is the time elapsed 
between successive peaks. It is evident that more regular 
the induced dynamics the lower the value of the com- 
puted NV. For purely periodic dynamics the NV goes to 
zero. 

Fig. Ufa) — (c) (left panel) show the time series of the 
floating potential fluctuations for different noise levels 
and Fig[7fd) (right panel) is the experimental NV curve 
as a function of noise amplitude D. The point (a) in 
Fig [Did) (time series shown in Fig. Ufa)) is associated 
with a low level of noise where the activation threshold is 
seldom crossed, generating a sparsely populated irregular 
spike sequence. As the noise amplitude is increased, the 
NV decreases, reaching a minimum (b) in Fig[7td) (time 
series shown in Fig. Efb)) corresponding to an optimum 
noise level where maximum regularity of the generated 
spike sequence is observed. As the amplitude of superim- 
posed noise is increased further, the observed regularity 
is destroyed manifested by an increase in the NV; label 
(c) in Fig[7Jd) (time series shown in Fig. 0(c)). This is a 
consequence of the dynamics being dominated by noise. 
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FIG. 4: Left panel: Timeseries showing that the period of relaxation oscillations increases with augmenting DV. Right panel: (a) 
Exponential increment of the time period (T) with DV and (b) In \V — Vh | vs T curve can be fitted by a straight line indicating 
an underlying homoclinic bifurcation. 
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FIG. 5: Evolution of glow size of the anode glow with increasing DV (a) - (h). 



V. DISCUSSION 

The effect of noise has been studied experimentally 
near a homoclinic bifurcation in glow discharge plasma 
system. Our study demonstrates the emergence of SR 
for periodic subthreshold square pulse signals and the in- 



duction of CR via purely stochastic fluctuations. In SR 
experiments, the efficiency of information transfer was 
quantified using AMD instead of the power norm which 
has been utilized elsewhere @. The advantage of using 
this method in comparison to the power norm (Co(0)) [5| 
lies in the fact that AMD remains independent of the lag 
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FIG. 6: Noise induced system response in the floating potential fluctuations for low, medium , and high amplitude noise in 
conjunction with a subthreshold periodic square pulse. The right panel [Fig [6] shows the AMD as a function of noise amplitude 
for the experiment performed at Vo — 307 V and pressure= 0.39 mbar. Left panel shows the subthreshold periodic pulse 
train and the three time series of floating potential fluctuations at low level noise (a); at optimum noise value (b) and at high 
amplitude noise (c). 
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FIG. 7: Emergence of coherence resonance: The right panel shows the NV as a function of noise amplitude for the experiments 
performed at Vo = 344 V and pressure = 0.5 mbar. Left panel: The time series of the floating potential fluctuations for (a) low 
noise (b) optimum noise and (c) high level noise 



between the measured floating potential and the applied 
periodic square pulse. This is of relevance in our exper- 
imental system, where invariably there exists a lag, at 
times varying in time due to the parameter drifts. For 
the CR experiments it was occasionally observed that 
while with an initial increase in noise amplitude (D) NV 
reaches a minimum, the subsequent rise of NV for even 
higher amplitudes of noise was suppressed. This leads 
to the modification of the unimodal profile, a signature 
of the CR phenomenon. A possible explanation for this 



suppression is that by virtue of the superimposed high 
frequency noise (bandwidth 500 KHz) and fast respond- 
ing internal plasma dynamics, the system has the capa- 
bility of exciting high frequency regular modes within the 
ion plasma frequency (105 kHz). This in turn leads to 
the persistence of low NV values. Finally, in Refs fl2l.lT3j 
both the destructive and constructive role of noise (CR 
only) have been reported for glow discharge and magne- 
tized rf discharge plasma systems respectively. However, 
both these experiments were carried out in the vicinity 
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of the Hopf bifurcation. In contrast, for the present work 
we studied both stochastic (SR) and coherence resonance 
(CR) in the neighborhood of the homoclinic bifurcation. 

APPENDIX A: CALCULATION OF ABSOLUTE 
MEAN DIFFERENCE(AMD) 

Regularity in the stochastic resonance has been done 
by calculating cross-correlation (C = | < \{xi — < x\ > 
)(x2— < x 2 >)] > |). But in this case this is not suit- 
able, because, as we also measuring floating potential at 
different location in side the plasm there is always a lag 
between periodic signal that is applied in the plasma and 
output. This lag also varies with time because the plasma 
conditions are changing continuously with time. There- 
fore cross-correlation between output and input signal 



gives wrong estimation. So we have proposed a statistics 
which will be independent of lag and defined as follows: 

1. First calculate mean inter-peak distance (S) of the 
periodic signal. 

2. Calculate inter-peak distances (t p ) of the output 
signal. 

3. Calculate {{t p - 5)/ 5) 

4. Take absolute. 

5. This will perfectly describe the regularity for 
stochatic resonance for periodic subthreshold signa. 

and this may be called Absolute Mean Difference 
(AMD). 

AMD = | < - 1) > | where t p , S arc the output 
signal inter-peak distance and mean peak distance of the 
subthreshold periodic signal respectively. 
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